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Amendments to the Claims 

This listing of claims will replace all prior versions, and listings, of claims in this application. 

The following amendments do not constitute an admission regarding the patentability of the 
amended subject matter and should not be so construed. Amendments to the claims were made for 
purposes of more clearly stating the claimed subject matter and do not add new matter or alter the 
scope of the claims. 

Listing of Claims: 

1. (Previously presented) A method for constructing a gene network of relationships 
between a plurality of genes of a set of genes and generating a graph representing the gene network 
of relationships between the genes comprising the steps of: 

(a) providing a quantitative time course data library for the set of genes of an organism, 
said library including expression results based on time course of expression of each gene in said set 
of genes, wherein the expression results comprise gene expression data at a plurality of time points 
for each gene in said set of genes, quantifying an average effect on expression of each gene at each 
time point by each other of said genes at each time point, and quantifying a measure of variability of 
the effect on expression of each gene at each time point on each other of said genes; 

(b) creating a gene expression matrix from said library, wherein the gene expression matrix 
comprises matrix data comprising the effects on expression of each gene at each time point by each 
other of said genes at each time point; 

(c) applying a Bayesian computational model to the matrix data, wherein said Bayesian 
model comprises minimizing a BNRCdynamk criterion; and 

(d) generating the gene network of relationships between said genes and generating a graph 
representing the gene network of relationships between said genes. 

2. (Canceled) 
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3. (Previously presented) The method of claim 1, wherein said step of minimizing a 
BNRCdynamic criterion comprises using a non-linear curve fitting method selected from the group 
consisting of polynomial bases, Fourier series, wavelet bases, regression spline bases and B-splines. 

4. (Previously presented) The method of claim 1, wherein said data library is created by 
obtaining the gene expression data at the plurality of time points from a time course study studying 
altering gene expression over time. 

5. (Previously presented) The method of claim 1, wherein said step of minimizing said 
BNRC^namic criterion further comprises selecting a Bayesian mode using a backfitting algorithm. 

6. (Previously presented) The method of claim 1, wherein said step of minimizing a 
BNRCdynamic criterion further comprises using Akaike's information criterion. 

7. (Previously presented) The method of claim 1, wherein said step of minimizing a 
BNRCdynamic criterion further comprises using maximum likelihood estimation. 

8. (Previously presented) The method of claim 1, wherein said plurality of genes are 
associated with a cell cycle of the organism. 

9. (Previously presented) The method of claim 1, wherein said measure of variability is 
variance. 

10. (Previously presented) The method of claim 1, wherein said non-linear curve fitting 
method is a non-parametric method. 

11. (Original) The method of claim 10, wherein said non-parametric method for 
minimizing a BNRCdynamic criterion includes using heterogeneous error variances. 

12. (Previously presented) The method of claim 11, wherein said step of minimizing a 
BNRCdynamic criterion further comprises the steps of: 

(1) making a score matrix whose (i, j) th element is the BNRCdynamic score of the graph 
gene,-* gene,; 
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(2) implementing one or more of adding the BNRC j ,fy„ nm , c score of the effect of parent 
gene, on gene, to a calculation of BNRQ^ M?m - c , removing the BNRC j dynamic score of the effect of 
parent gene, on gene,- to a calculation of BNRC^^jc, and reversing the BNRC j (fy „ am , c score of the 
effect of parent gene ; on gene/ to a calculation of BNRC rf> ,„ flm;Cj which provides the smallest 
BNRCdynamic', and 

(3) repeating step 2 until the BNRC^^/c does not reduce further. 

1 3 . (Original) The method of claim 1 1 , wherein said step of minimizing a BNRCrfy„ amiC 
criterion further comprises the step of applying a hill-climbing algorithm to minimize BNRd dynamic . 

1 4. (Previously presented) The method of claim 1 1 , wherein an intensity of an edge is 
determined using a bootstrap method. 

15. (Previously presented) The method of claim 14, wherein said bootstrap method 
comprises the steps of: 

( 1 ) providing a bootstrap gene expression matrix by randomly sampling a number of 
times, with replacement, from the original gene library expression data; 

(2) estimating the gene network for genej and genej; 

(3) repeating steps ( 1 ) and (2) T times, thereby producing T gene networks; and 

(4) calculating the bootstrap edge intensity between gene; and genej as (ti+t 2 )/T. 
16-17. (Canceled) 

1 8. (Previously presented) The method of claim 1, comprising grouping said genes into 
one or more equivalence sets. 

19. (Canceled) 

20. (Previously presented) A method for constructing a gene network model of a system 
containing a network of genes from time course gene expression data, said method comprising: 
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applying a Bayesian computational model to the time course gene expression data, 
wherein applying said Bayesian computational model comprises minimizing a BNRCdynamic 
criterion; and 

generating a graph representing the gene network model resulting from applying the 
Bayesian computational model to the time course gene expression data wherein BNRC dynamic 
criterion is minimized. 

2 1 . (Original) The method of claim 20, wherein minimizing the BNRCdynamic criterion 
comprises using a non-linear curve fitting method selected from the group consisting of polynomial 
bases, Fourier series, wavelet bases, regression spline bases and B-splines. 

22. (Previously presented) The method of claim 20, wherein minimizing the BNRCdynamic 
criterion comprises selecting a Bayesian model using a backfitting algorithm. 

23. (Original) The method of claim 20, wherein minimizing the BNRCdynamic criterion 
comprises using Akaike's information criterion. 

24. (Original) The method of claim 20, wherein minimizing the BNRCdynamic criterion 
comprises using maximum likelihood estimation. 

25. (Original) The method of claim 20, wherein minimizing the BNRCdynamic criterion 
comprises using a non-linear curve fitting method, wherein the non-linear curve fitting method is a 
non-parametric method. 

26. (Original) The method of claim 25, wherein the non-parametric method includes 
using heterogeneous error variances. 

27. (Previously presented) The method of claim 26, wherein minimizing the BNRC^fl^c 
criterion further comprises the steps of: (1) making a score matrix whose (i, j) th element is the 
BNRCdynamic score of the graph gene, -> gene,; (2) implementing one or more of adding the 
BNRCdynamic score of the effect of parent gene, on gene y - to a calculation of BNRCdynamic, removing 
the BNRC* dynamic score of the effect of parent gene, on gene/ to a calculation of BNRCdynamic, and 
reversing the BNRC } dynamic score of the effect of parent gene,- on gene, to a calculation of 
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BNRCdynamic which provides the smallest BNRdynamic', and (3) repeating step 2 until the 
BNRCdynamic does not reduce further. 

28. (Previously presented) The method of claim 26, wherein minimizing the BNRCdy namic 
criterion further comprises the step of applying a hill-climbing algorithm to minimize the 

BNRC 0) dynamic . 

29. (Previously presented) The method of claim 26, wherein an intensity of an edge is 
determined using a bootstrap method. 

30. (Previously presented) The method of claim 29, wherein said bootstrap method 
comprises the steps of: 

( 1 ) providing a bootstrap gene expression matrix by randomly sampling a number of 
times, with replacement, from the original gene library expression data; 

(2) estimating the gene network for genej and genej; 

(3) repeating steps (1) and (2) r times, thereby producing Tgene networks; and 

(4) calculating the bootstrap edge intensity between genej and genej as (ti + t 2 )/T. 
31-44. (Canceled) 
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